Optical conductivity for a dimer in the Dynamic Hubbard model 



G. H. Bach and F. Marsiglio 
Department of Physics, University of Alberta, Edmonton, Alberta, Canada, T6G 2G7 

(Dated: March 2, 2013) 

The Dynamic Hubbard Model represents the physics of a multi-band Hubbard model by using 
a pseudo-spin degree of freedom to dynamically modify the on-site Coulomb interaction. Here 
we use a dimer system to obtain analytical results for this model. The spectral function and the 
optical conductivity are calculated analytically for any number of electrons, and the distribution of 
optical spectral weight is analyzed in great detail. The impact of polaron-like effects due to overlaps 
between pseudo-spin states on the optical spectral weight distribution is derived analytically. Our 
conclusions support results obtained previously with different models and techniques: holes are less 
mobile than electrons. 



I. INTRODUCTION 

The occurrence of electron-hole asymmetry in tunnel- 
ing spectra^ and the anomalous behaviour in the optical 
conductivity sum rule at the superconducting transition 
temperature^^ both contribute to the possibility that 
the superconductivity in the cuprate materials is unusual 
in several respects. In particular, the notion of 'kinetic 
energy-driven superconductivity' is implied by the opti- 
cal experiments, as predicted almost ten years in advance 
of these experiments.^ 

These experiments indicated that a significant trans- 
fer of spectral weight occurs in the cupratesP both in 
the normal state as a function of temperature, and as 
a result of the superconducting transition. More impor- 
tantly, perhaps, is the range of frequencies affected by 
the transition, as a significant amount of spectral weight 
is transferred from very high frequencies to very low fre- 
quencies. Thus, there is an apparent violation of the 
Ferrell-Glover-Tinkham optical sum rule, as an examina- 
tion of the low frequency region alone shows a spectral 
weight discrepancy. This indicates that physics beyond 
the usual paradigm of "energy lowering due to poten- 
tial energy considerations" is at work; in particularly the 
anomalous sign of the change in low frequency optical 
spectral weight indicates that some mechanism involving 
the kinetic energy of the carriers is at work. 

Earlier modelling of cuprate superconductivity^ in- 
cludes some of this physics — this is what motivated 
the relatively early theoretical discussion of optical spec- 
tral weight transfer — but a more recent theoretical 
model, advanced more than ten years ag goes further 
to explain some of the anomalous features in the spec- 
troscopic measurements of the cuprates; this is the so 
called Dynamic Hubbard model. This model utilizes a 
phenomenological pseudo-spin degree of freedom at each 
lattice site designed to mimic orbital relaxation effects 
which necessarily occur in real atoms. As far as opti- 
cal spectral weight transfer is concerned, this model in- 
cludes higher frequency excitations (here modelled by the 
pseudo-spins), and therefore, while we do not address su- 
perconductivity or temperature effects in this paper, by 
using the Dynamical Hubbard model, we can study how 
spectral weight transfer occurs as a function of doping. 
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where the pseudo-spin degree of freedom is here repre- 
sented by a Pauli operator Oi at each site; it interacts 
with the electron charge through the double occupancy 
term, and contributes a dynamical interaction in addi- 
tion to the usual Hubbard interaction. The rest of the 
Hamiltonian is as follows: the first term represents the 
electron kinetic energy within a tight-binding model with 
one orbital per site. Note that we are really trying to un- 
derstand physics that originates in processes involving 
multiple orbitals. It is desirable to minimize the com- 
plexity by retaining a single orbital, and it is then the 
pseudo-spin that acts to mimic the physics of carriers 
undergoing transitions between multiple orbitals when 
the local occupation changes. The second term determines 
the electron density through the chemical potential, the 
third term defines the two level system for the pseudo- 
spin degree of freedom at each site, and the last term 
is the on-site interaction which, in addition to the short 
range Coulomb repulsion represented by U, is also mod- 
ulated through a coupling constant gujQ by the state of 
the pseudo-spin. 

When the double occupancy is high, the pseudo-spin 
will reside in its excited state for the sake of minimiz- 
ing the Coulomb repulsion, much like the phenomenon 
in real atoms, where two electrons will sacrifice having 
a minimal electron-ion energy and spread out amongst 
the excited orbitals in order to minimize their Coulomb 
repulsion. In the opposite limit, when the double oc- 
cupancy is very low, electrons will tend to stay in the 
lowest energy state available in the given atom (loosely, 
the Wannier state which is being modelled in the tight- 
binding Hamiltonian), and the Coulomb energy will be 
high, though irrelevant, since only rarely will two elec- 
trons occupy the same site. 

The Dynamic Hubbard model contains at leas t som e 
of the phenomenology of hole superconductivity, 121 ^^^ 
proposed more than twenty years ago. In particular, 
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the model contains electron-hole asymmetry, where holes 
at the top of a band are heavier than electrons at the 
bottom of a band. Electron-hole asymmetry can arise 
in a number of ways: just having further than nearest 
neighbour hopping can result in a band mass asymmetry, 
upon which polaronic mass asymmetry can buildP^ Even 
with just nearest neighbour hopping, lattice geometry 
can also result in electron-hole asymmetry, and again be 
amplified through polaron effects.^ Here, the electron- 
hole asymmetry arises through realistic interactions, and 
is strongly connected to the fact that holes pair more 
readily than electrons by lowering their kinetic energy in 
the superconducting phase, a phenomenon supported by 
the anomalous observations in the optical conductivity 
sum rule mentioned above. More recently, the Dynamic 
Hubbard model has been explored with Dynamical Mean 
Field Theory (DMFT) P 1 In particular, the spectral func- 
tion and the optical conductivity were calculated, in the 
normal state, to illustrate the electron-hole asymmetry 
present in the model. 

In this paper we will focus on the optical conductiv- 
ity, and provide a complementary calculation involving a 
simple dimer. Such a small system does not constitute a 
very realistic system; however, dimer calculations have an 
illustrious history for providing insight into models E2HHD 
Furthermore, following Ref. [7], a dimer calculation pro- 
vides good insight into the processes that contribute to 
the conductivity. We will briefly review the optical con- 
ductivity and the sum rule in the next section, and de- 
scribe the details of the dimer calculation in the third 
section, in perturbation theory; we have also performed 
exact diagonalizations to delineate the regime of validity 
of the perturbative results. This is followed with a dis- 
cussion of the results, particularly in light of the DMFT 
results reported in Ref. [TH]. We then conclude with a 
summary. 



II. OPTICAL CONDUCTIVITY AND RELATED 
SUM RULES 

The real part of the optical conductivity cr 1 (w), as a 
function of frequency, u>, at finite temperature T can be 
written as: 
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where |m) and \n) are eigenstates of the Hamiltonian of 
the system, Z is the partition function, /3 = 1/ksT and J 
is the current operator, obtained through the polarization 
operator PHI As formulated by Kubo in 1957,^ the optical 
conductivity satisfies the general sum rule, 
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where r denotes the type of charge carrier, and n r , e r 
and m r are the number density, charge, and mass, re- 
spectively, of the r-type carrier, and /i, v are the indices 
of the conductivity tensor. For an isotropic electron sys- 
tem, the sum rule (|| is rewritten as: 
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with m the bare mass of electrons and n the total electron 
density.^ 

In condensed matter systems we often work with ef- 
fective Hamiltonians, for example formulated for a single 
band within tight-binding. One can then formulate a 
sum rule restricted to that single band, and obtairP^U 
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where is the single electron occupation number, and 
£fc is the dispersion relation for the non-interacting elec- 
trons. In reality, the integration in Eq. ([5]) is taken up 
to a cut-off frequency uj c determined experimentally in 
order to allow only intraband transitions, and to avoid 
the inevitable interband transitions which are not part of 
the sum rule Eq. Theoretically, the right-hand-side 
(RHS) of Eq. |5| is often usedpSlas this is much simpler 
to calculate. Furthermore, when only nearest neighbour 
hopping is allowed on a hypercubic lattice, the sum rule 
Eq. (p3| reduces toPl 
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with K = —tJ2(ij) a= j ti( c lcr c ja- + h.c.) and a the lattice 
constant. From the RHS of Eq. ffify, it is clear that the 
optical sum depends not only on the external parame- 
ters (such as the temperature) but also on the electronic 
structure of the system. The validity of using the kinetic 
energy on the RHS instead of the expression in Eq. (|5j, 
even when the dispersion is not just nearest neighbour 
hopping has been explored in Refs. [28l and [27] . to which 
the reader is referred. In the following we assume Eq. © 
holds. 

Measurements of the optical conductivity sum rule in 
a number of the cuprate superconductors generally show 
an increase of spectral weight in the low frequency regime 
in the superconducting state, at least in the underdoped 
and optimally doped materialsP^ This enhancement of 
the optical sum at the superconducting transition tem- 
perature conflicts with the result from BCS-like super- 
conductivity where an increase in the kinetic energy (and 
therefore a decrease in the single band optical sum) is ex- 
pected instead. This can be explained through a number 
of different scenarios, examples of which are preformed 
pair^H, and phase fluctuations pSHU I n contrast to the 
model considered here, many of these calculations have, 
as a key in gred ient, proximity to a near by M ott insu- 
lating state!^J22l Note that some authors^ 1 ^ 1 have at- 
tributed the anomalous temperature dependence of the 
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low frequency optical spectral weight to a cutoff effect (re- 
quired in the experimental analysis). Karakozov et alP^ 
have also attributed the anomalous change in spectral 
weight at the superconducting transition temperature to 
a cutoff effect, though this has been refuted in Ref. [25] . 
See Ref. [35] for a brief review. 

On the other hand, models like the hole mechanism 
of "kinetic energy driven" superconductivity support the 
idea of minimizing the total energy by reducing the ki- 
netic energy, and therefore the optical sum has an anoma- 
lous temperature dependence below T c Pf When the sys- 
tem goes superconducting, the missing optical spectral 
weight is predicted to be distributed over the whole range 
of frequencies, i.e. weight is transferred from high fre- 
quency to low frequency and the lo w en ergy sum rule 
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appears to increase as a consequence. 1 



III. OPTICAL CONDUCTIVITY IN A DIMER 

We proceed now with a brief discussion of the site 
Hamiltonian followed by a detailed description of the 
dimer. 



A. Properties of the Hamiltonian 

Following Refs. [TBI and I3"8] . we begin with the on-site 
Hamiltonian for electrons: 
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dhm = uoo-'x + guj o- l z + [U - 2guj a l z }n lt n zl . (7) 



The solutions are provided in detail in Refs. [TBI and 138] : 
for n electrons the ground state (|n)) and the first excited 
state (|n)) are 



with 
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The eigenvalues (ground state e(n) and excited state 
e(n)) are 



e(n) = S ny2 U - u) Q \Jl + g 2 



Especially important for the hopping processes is the 
overlap of background spin states with different numbers 
of electron; these are 
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These parameters play an important role for the spec- 
tral function; these are defined in Ref. [35] as A n+l n for 
electron destruction in a system of n + 1 electrons (and 
A n ^ n+ i for electron creation in a system of n electrons. 
For a single site these single particle spectral functions 
are^l 

Aio(w) =A i(u) = 6(u) 

A 12 (u) = A 21 (w) - S 2 S{uj) + (1 - S 2 )5(lj - fio).(13) 

Even though the spectral weight is calculated for a single 
site, it is clear that there is a reduction of the weight 
at zero frequency if the second electron is added to 
the one-electron ground state. The reason is because 
there are two possibilities: the pseudo-spin can remain 
in the same state as the first electron with a probability 
S 2 < 1, or it can become excited with an energy cost 
fio = 2wo yl + g 2 ■ In the thermodynamic limit, this ef- 
fect is known as the reduction of quasiparticle weight by 
transferring part of the coherent contribution (at w = 0) 
to the incoherent part (at large ui), resulting in a one 
particle spectral weight, z < 1. Since the quasiparti- 
cle weight is inversely proportional to the effective mass, 
z ~ m/m*, this statement means that holes are heavier 
(or more 'dressed®') than electrons. 

For calculating the optical conductivity in perturba- 
tion theory, the Hamiltonian is divided into two parts: 

H = H Q + H', (14) 

H Q = ^2(u} of + guoo-? ) + (U- 2guj o-?)n it riii ( 15 ) 

i 

H' = K = -t ]T(4c^+ Cj W), (16) 

where H is the site Hamiltonian and W is the hopping 
part which is considered as a perturbation under the fol- 
lowing conditions. Based on the definition of the over- 



laps between the pseudo-spin ground states in Eq. (12 1, 



we can define a pseudo-spin state for a given number of 
electrons in terms of the eigenstates involving a different 
number of electrons. These overlaps contain the back- 
ground deformations (modelled by the pseudo-spin) that 
must be 'dragged' along as the electron hops. Thus, fol- 
lowing Ref. [7]: 
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where 



S = \j\ - S 2 = 



(21) 



We wish to solve this problem in all number sectors (one, 
two, and three electrons). We cover these in the following 
subsections. 



B. Three electron sector 

Beginning with three electrons on the dimer (a hole- 
like configuration), the (non- normalized) ground state 
wave function is given in first order perturbation theory 
as: 
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are the 8 basis states required to span the Hilbert space 
for the three electron sector. Here the subscript 'e' refers 
to the even states, and a and b are the indices of the 
first and second site, respectively, in the dimer. Note 
that kets followed by the subscript a or b have numbers 
0, 1, or 2 (with or without bars on top) that refer to the 
pseudo-spin eigenstates defined in Eq. pi), whereas kets 
followed by e (for 'even') or o (for 'odd 7 ) refer to linear 
combinations of product states of electrons and pseudo- 
spin eigenstates, as, for example, in Eq. (23 1. 



The three-particle ground state energy, to first order 
in the hopping perturbation, is given by: 
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and the excited state energies for the three electron sec- 
tor, are, 
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FIG. 1. Schematic depiction of optical transitions in a dimer 
with three electrons. The lines (both solid and dashed) show 
the two available levels of the pseudo-spin energy at each 
site; the solid lines correspond to the occupied pseudo-spin 
state and the dashed lines correspond to the unoccupied state. 
Transitions between states with the same pseudo-spin energy 
levels are diagonal; these contribute to the intraband con- 
ductivity, while non-diagonal transitions between states with 
different pseudo-spin energy levels modify the interband con- 
ductivity. State labels are those found in Eq. (231 in the text, 
where they are given in full even or odd form. 



Note that there are degeneracies at zeroth order between 
even and odd states; these are broken in first order per- 
turbation theory, as is explicitly written in Eq. ( 24 ) and 



(23) the first of Eqs. ( 25 1 . Only the zeroth order energies (and 



wave functions) are needed for the other excited states, 
and that is what is written here. Also, the Hamiltonian 
will only couple states of a given parity, whereas the con- 
ductivity will couple only states of opposite parity. 

The optical conductivity for the dimer at zero tem- 
perature can be calculated for the three electron sector 
as: 
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where |m) e are the excited states of the system (only even 
parity is required since the ground state has odd parity 
— these are given in Eq. (23)), and 
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is the current operator. By acting with the J operator on 
the ground state IV^)' we connect to three of the excited 
states (all even parity) given in Eq. ( p3| . Note that since 
the current operator is already of order t, only zeroth 
order wave functions are required, but for the first even 
parity excited state, first order corrections to the energy 
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are required in the denominator of Eq. ( 26 ) to break the 
degeneracy. 

Operating with the current operator on the unper- 
turbed ground state gives 

J|l)o = l -f (s 2 \l) e - V2SS\3) e - S 2 |4) e ) (28) 

so that the optical conductivity for three electrons in- 
cludes three peaks. These three transitions are shown 
schematically in Fig. 1; the analytical expression for the 
optical conductivity (for three electrons) is 
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The optical sum rule can be checked by calculating the 
expectation value of the K operator in the ground state, 
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which is precisely the combination of weights given in 
Eq. (29 1. The first contribution comes from intraband 
transitions — this would correspond to the Drude weight 
for an extended system. This Drude response is, however, 
weighted by the overlap S = (1|2) between the respec- 
tive ground states of the pseudo-spin with one and two 
electrons. This is referred to as a 'diagonal' transition 
in Fig. 1, since the background (here, the pseudo-spin) 
doesn't become excited in the transition. The second and 
the third peaks involve transitions corresponding to one 
and two pseudo-spin excitations, f^o; these are recognized 
as interband transitions in the language of multiple band 
models. This first order perturbation approximation re- 
sult remains valid as long as t/Q <C 1 Comparisons with 
exact results will be shown below. 



C. Two electron sector 

The same procedure can be performed with the more 
difficult case of two electrons. In this case there are 16 
basis states, and again they can be divided into 8 even 
and 8 odd states. We use a slightly different notation 
- there are now states involving double occupation of 
a single site, and those involving only single occupation. 
The states with double occupation are 
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where the 'd' in the ket stands for 'double', and the 
subscripts '1,2,3,4' simply enumerate the states, start- 
ing with the lowest, |di), where both pseudo-spins are in 
their ground states, The three other basis states corre- 
spond to the first site have an excited pseudo-spin state, 
the second site having an excited pseudo-spin state, and 
both sites having excited pseudo-spin states, respectively. 

Similarly, for the basis states involving only singly oc- 
cupied sites, the 's' in the ket stands for 'single'; these 
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The first has both pseudo-spin states in the ground state, 
with the three others having excited pseudo-spin states 
as in the case of the doubly occupied states. Note that an 
equally valid set of states combines \s2) and |ss) symmet- 
rically and anti-symmetrically (as we did in the middle 
two basis states of Eq. ( 23 ) for the three electron sector) . 

We want the unperturbed ground state to reside in the 
space of states involving only ground state pseudo-spin 
states. Confining ourselves to this sector only, the ground 
state wave function is 
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This will be an accurate ground state wave function as 
long as the pseudo-spin excitation energy remains the 
largest energy scale in the problem, i.e. U <§; £!o, along 
with the restriction already used, t <C Sl - When U <C 
2tS the ground state consists of nearly equal amplitudes 
of the two basis states; on the other hand, when U 3> 2tS, 
ao ~ and bo ~ 1, and the singly occupied basis state 
dominates the ground state, as expected. With these 
same assumptions the two electron ground state energy 
is given as 



4 2) « -flo + U/2 - ^/(U/2) 2 + (2tS) 2 . 



(36) 



To determine the optical conductivity we need the re- 
sult of the current operator on each component of the 
ground state; the result is 
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FIG. 3. The normalized optical conductivity as a function of 
frequency in a dimer obtained by perturbation theory with 
U = 0, g = 3 and ujq = 1. The 5— functions in the formulas in 
the text are represented here as broadened Lorentzians (with 
a width 0.05). Note the decrease in relative low frequency 
spectral weight (labelled as pi) as one goes from one to two to 
three electrons, indicating a reduced mobility as the number 
of electrons increases. Also note as the number of electrons 
decreases the higher frequency relative weight decreases: for 
the two electron case spectral weight is entirely absent at 2Qo 
(labelled as ps) while for the one electron case weight is absent 
even for cj = flo (labelled as P2)- 



FIG. 2. Schematic depiction of the optical conductivity with 
two electrons in a dimer. Parts (a) and (b) refer to the two 
basis states that make up the ground state wave function 
given by Eq. ( |33[ ) with ground state energy given by Eq. ( |36[ ). 
The first has only (two) non-diagonal transitions (the second 
transition, to states given by |s2)e, is qualitatively the same as 
the one shown), while the second basis state has both diagonal 
and non-diagonal transitions, as shown. The state represen- 
tations are schematic only; formulas in the text represent the 
full (even or odd) state. 



Fig. 2 summarizes these transitions schematically. The 
unperturbed energies associated with these states are 
readily determined by inspection from Eqs. (31 ) and (32 1. 



Using the analogue of Eq. ( 26 1 appropriate to two elec- 



trons, with the ground state now a linear superposition 



of basis states given by Eq. ( 33 1 , the optical conductivity 
for two electrons is obtaincc 
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(39) 
0) or 



plays the role of a low energy scale, i.e. t (as U 
U (as * -> 0). 

The two electron optical conductivity also has three 
peaks, as seen in Eq. (381 or schematically in Fig. 2. In 



the perturbative approach we have used, however, the 
two higher frequency peaks should be viewed as one oc- 
curring at a frequency scale of order Hq that has been 
split by a low energy scale of order U or t. The analog of 
the third peak in the 3 electron case is absent here; there 
is no peak at 2f2o- 

The first (low energy) peak corresponds to the diagonal 
transition from state |si) to \di) , where, as in the three 
electron case, the hopping of an electron between sites oc- 
curs without modification of the pseudo-spin background. 
This diagonal transition corresponds to the Drude-like 
(or coherent) part of the optical conductivity, though it 
may extend to a range of (low) frequencies on the scale 
of U as well as t (as represented by e). 

The second peak at frequency lu = Q,q — U + e is given 
by the two transitions: \di) e to both |sa)e and |s3) e . The 
third peak at frequency f2 + e is obtained by the \si) a to 
|d3) state transition, as detailed in Fig. 2. The second 
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FIG. 4. The dependence of the various normalized spectral weight contributions to the optical conductivity on electron-pseudo- 
spin coupling g for three electrons (a) and two electrons (b). In (c) and (d) we show the same quantities as a function of wo, 
while in (e) and (f) they are shown as a function of U. Perturbation results are shown as curves, while exact results (for the 
dimer) are shown by symbols, as indicated. Note that usually the Drude weight dominates; however, for sufficiently large g 
(and U in the case of two electrons) the Drude weight is significantly reduced, indicative of reduced mobility, especially for the 
highly (electron) doped regime. 



and third peaks correspond to non-diagonal transitions, 
which means that the pseudo-spin background is excited 
in the transition; this in turn corresponds to transitions 
involving higher energy bands, not included in our start- 
ing Hamiltonian. There is also an overall factor of two 
enhancement because there are now two carriers instead 



of the one carrier present in the three electron case pre- 
sented above (or the one electron case shown below). 



One can again verify the conductivity sum rule; for 
this we need the ground state wave function given to 



first order in t/Clo- Straightforward calculation^ gives 

(2) 

\% Oo|dl) e + &o|si)o + 

_^ (|8s)0 + w 

and an evaluation of the kinetic energy expectation value 
gives 

-(^\K\4 2) ) = 

V(J7/2) 2 + (2tS) 2 ^o-U + e O + e ' 



again in agreement with the weights in Eq. ( 38 ) 



D. One electron sector 

Finally, similar calculations in the one electron sector 
are particularly simple, because the single site pseudo- 
spin ground states for zero and one electron are identi- 
cal. This means that there are no pseudo-spin excitations 
arising from application of the current operator to the 
single electron ground state. The optical conductivity 
for one electron is given by the simple expression, 



ne/t 



8{w-2t/h), 



(42) 



which contains only a Drude contribution, with no nor- 
malization (as required by the sum rule). 



IV. DISCUSSION 

We wish to show the relative contributions to the con- 
ductivity for the three different electron densities that we 
can access. Because this is a tight binding model it will 
not conserve total oscillator strength. This makes com- 
parisons for different numbers of electrons and/or differ- 
ent parameter values difficult!^ Here, for a given number 
of electrons, we will normalize the conductivity to the 
overall spectral weight in the conductivity for that num- 
ber of electrons. Fig.[3]shows the two-site optical conduc- 
tivity with one, two and three electrons, using U = and 
g = 3, with u! = 1, as a 'standard' set of parameters. 
While not necessarily realistic, they are chosen specifi- 
cally to remove Mott complications at half filling; here 
the presence of U will result in a significant decrease in 
low frequency spectral weight, in spite of the increase in 
number of available carriers. These parameters will serve 
to illustrate the spectral weight transfer physics inherent 
in this model. At the same time, it is clear that the per- 
turbation calculation is valid if loq is large enough so that 
Oo is much larger than U, t as assumed above. 

The simplest case is clearly that of one electron. 
Eq. ( 42 ) and Fig.[3]show that a single peak is present; it is 
located at to = 2t/H, a non-zero value only because we are 
dealing with a dimer, and not the thermodynamic limit. 



In the thermodynamic limit this would be a ^-function 
at zero frequency, representing the Drude contribution. 
Normally this peak would be broadened through, for ex- 
ample, impurity scattering, but here (and even in the 
DMFT study of Ref. [T^]) it remains a (^-function, broad- 
ened in Fig. 3 artificially by hand, so as to be visible. For 
small but non-zero electron densities U would contribute 
as well, but for the most part this picture would remain 
unchanged. In particular, excited pseudo-spin states are 
essentially absent (see Fig. 3 in Ref. [16], which shows 
that that the expectation value of the pseudo-spin op- 
erator is essentially its ground state value for low den- 
sities). Because the electron density is relatively dilute, 
there are very few optical transitions involving doubly oc- 
cupied sites, and therefore it is not possible to excite the 
pseudo-spin excited state. For the dimer the only repre- 
sentative sector for this physics is the one electron sector 
(two electrons already constitutes a somewhat crowded 
lattice) . 

In contrast the three electron sector represents the 
most 'crowded' situation for a dimer, while the two elec- 
tron sector is somewhat in between, and, as mentioned 
previously, the absence of a Coulomb repulsion (U = 0) 
aids to highlight the pseudo-spin physics, and suppress 
the Mott-related physics (which, from our point of view, 
is not essential, and will complicate the analysis). Refer- 
ring to Fig. 3, note that the 3 electron conductivity has a 
significant relative contribution at high frequency (20o)j 
this is entirely absent in the two electron conductivity — 
it has been pushed down to lower frequency (Oq). The 
reason for this is as follows: with three electrons, the 
ground state consists of a doubly-occupied and a singly- 
occupied site, with the respective pseudo-spins at each 
site in its ground state — see the first of Eqs. (23 1. An 



optical transition can result in one of the three states 
shown in Fig. 1; one of these, state |4) e , has two ex- 
cited pseudo-spin states, corresponding to an energy 2il - 
One of these excitations comes from the site with a sin- 
gle electron — before the transition this site was doubly 
occupied, and the ground state for this configuration re- 
quired a pseudo-spin ground state corresponding to two 
electrons. Since one has left, there is now a component 
of the pseudo-spin which corresponds to an excited state 
for the one electron configuration. Similar remarks ap- 
ply for the site that was previously singly occupied and 
is now doubly occupied. For two electrons this cannot 
happen — see Fig. 2 and note the absence of an alter- 
native involving two excited pseudo-spin state. This is 
because the pseudo-spin ground state is the same for an 
empty and singly-occupied site — see Eq. |9]) or the first 
of Eqs. (pi where T = 1. 



This accounts for the peak structure for the various 
electron sectors in Fig. 3. The dimer calculations have 
an 'all or nothing' character to them — no high frequency 
(2fio) peak for the two electron sector, and not even an 
intermediate (f2o) peak for the one electron sector. Of 
course in the DMFT calculations the changes from one 
electron density to another are smoothed out, as one can 
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see in Fig. 13 of Ref. [16 . The other feature that is 
apparent in Fig. 3 is the decrease of spectral weight in 
the relative Drude (low frequency) portion as one goes 
from the one electron to the two and then three electron 
sector. This is due to the polaron-like hopping renor- 
malization already discussed. The relative weight of the 
Drude portion is indicative of the coherence of the carri- 
ers, so again, in the dilute limit, electrons can hop while 
the background pseudo-spin degree of freedom remains 
in the same ground state at both the site from which the 
electron hops, and at the site to which the electron hops, 
because only empty or singly occupied sites are involved. 
In the more crowded lattice limit (here represented by 
the three electron sector), doubly occupied sites are nec- 
essarily involved, and then the background pseudo-spin 
has to adjust according to whether a singly or doubly 
occupied site is involved. 

The progression of spectral weight with electron num- 
ber is perhaps best exemplified by examining the conduc- 
tivity formulas, Eqs. ( 29|38|42 | for three, two, and one 
electron(s), respectively, for U — 0. Then the low fre- 
quency spectral weights are (omitting the common factor 
Tre 2 t/(2h 2 )) S 2 for three electrons (but one hole carrier), 
2S/2 = S for two electrons, and unity for one electron; 
these weights steadily increase by reducing the number 
of electrons, since S < 1 always, and this illustrates the 
principle that holes are less mobile than electrons. 

We analyze in more detail our results for the frequency 
dependence of the optical conductivity. The three elec- 
tron optical conductivity has three distinct peaks from 
low to high frequency, one at uj w t, one at uj « fioj and 
one at uj ~ 20o ; whose weights we denote pi, P2 and P3, 
respectively. In the two electron conductivity, there are 
again three peaks, but as explained above, the two high 
frequency ones are at the same characteristic frequency 
(identical if U = 0), so we will combine the weight from 
these two and denote it as P2; we will continue to use 
Pi for the lowest frequency peak, and of course for the 
one electron conductivity, there is only a low frequency 
Drude- like peak, which we will also denote as pi. In 
Fig. [4] we plot these weights to show how the spectral 
distribution of the optical conductivity varies with the 
strength of coupling g for 3 (a) and for 2 (b) electrons; 
in (c) and (d) we show the corresponding results as a 
function of uj , and in (e) and (f) results are shown for 
a variation of U. In all cases, the optical conductivity 
has been normalized to the total spectral weight for the 
parameters used in Fig. |3j separately for each electron 
number. 

As expected, increasing the coupling strength g be- 
tween the electron and the background (pseudo-spins) 
reduces the mobility of the electron as obtained in the 
spectral weight of the first peak p! (which would corre- 
spond to the Drude weight for an infinite lattice). This is 
simply due to the polaron effect mentioned above; with 
increased coupling, the amount of 'background adjust- 
ment' required as the electron hops increases. Physi- 
cally, the actual coupling in a given lattice is given by 



the amount of multi-orbital involvement required to min- 
imize the energy locally when two electrons try to accom- 
modate one another on the same site. Since we model this 
process with the pseudo-spin degree of freedom, we span 
a considerable parameter range in the figures. The abso- 
lute weight P2 of the second peak for the three electron 
case is given analytically as 4 ^ 1+ ^ 2 ^ ^ (see Eq. (29 1, 

and achieves its maximum value at g = \/2/3 w 0.8, 
which is independent of U. For the two electron case, 
the off-diagonal transition contributions, represented by 
P2, are quite negligible compared with the low frequency 
weight pi at weak coupling, but they play a more im- 
portant role at strong coupling. Note that results arising 
from a complete diagonalization of all the dimer states 
are also shown, and, for these parameter regimes, the 
agreement is excellent, as expected. The transitions de- 
noted by p2 and P3 represent incoherent processes; they 
may well correspond to the mid-infrared band that seems 
to feature so prominently in a wide variety of cuprate 
superconductors!^ The other experimental feature to 
which we can make contact with these dimer calcula- 
tions is the dependency on doping. Experimentally, the 
anomalies at the superconducting transition are most 
pronounced in the low hole regimeP consistent with the 
fact that the pseudo-spin physics in these calculations 
plays a large role precisely in this regime as well. 

Comparison with the results obtained from DMFT 
calculations^ is also possible. For example, in Fig. 13 
of Ref. [16], we show the conductivity as a function of 
frequency for various electron densities. Note that the 
parameters used in the DMFT calculation are in the more 
weak to intermediate coupling regime. Nonetheless, the 
calculations here are semi-quantitatively consistent with 
those. The first panel there refers to the very dilute limit 
(n = 0.1), and as suggested here, there is a single low fre- 
quency Drude peak. Of course there it is centred around 
zero frequency, while here it is at 2t, for reasons already 
explained. In the last panel in the same Fig. 13, n = 1.9, 
corresponding more to our present three electron calcu- 
lation. We expect weight at Qq = 2uj y/l + g 2 5.7, 
which is very close to the one shown there. Furthermore 
the Drude-like peak has reduced spectral weight (clearer 
in Fig. 15 of Ref. [TB]) compared to the result at n = 0.1. 
The expected peak at 2£lo is, however, barely present, 
and at a higher frequency than expected. It is not clear 
what the cause of this latter discrepancy is, especially in 
light of the quantitative accuracy of the other peaks. 

For completeness we have included plots to show the 
variation with ujq and U, where the expected behaviour 
occurs. Note that at half filling (two electrons) the exact 
results differ considerably from the perturbation theory 
results, as Mott physics becomes more prevalent (this is 
not surprising since this was not considered in the pertur- 
bative approach we took) . As ujq increases the results for 
three and two electrons become dominated by the Drude- 
like peak near the origin. Again, this is entirely expected, 
since pseudo-spin excitations become more and more en- 
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ergetically costly, and so, as seen explicitly in our per- 
turbative expressions, energy denominators increasingly 
suppress these transitions requiring excited pseudo-spin 
states, so that these play much less of a role as wo in- 
creases. As is clear from panels (c) and (d), exact di- 
agonalization results support these perturbative calcula- 
tions. 



V. CONCLUSIONS 

We have investigated spectral properties of the Dy- 
namic Hubbard Model on a dimer, primarily to gain 
a qualitative understanding of the physics of electron- 
hole asymmetry, and polaron-like mobility inherent in 
real atoms. Primarily we have investigated the spectral 
features of the optical conductivity with different num- 
bers of electrons. The physics we are trying to capture 
is that when electron movement results in a change from 
a doubly occupied site to a singly occupied site, or vice- 
versa, a considerable amount of 'background' adjustment 
needs to take place. In real atoms this is apparent in 
that the orbitals occupied by a single electron are con- 
siderably modified when two electrons occupy that same 
orbital. In the Dynamic Hubbard model, these modifica- 
tions are simulated by a pseudo-spin degree of freedom, 
at each site; an excited pseudo-spin state corresponds to 
an electron (partially) occupying an orbital that does not 
minimize the electron-ion energy, but does minimize the 
(local) electron-electron repulsion. 

Such processes will impact the optical sum rule; in 
particular, weight will be transferred over a considerable 
range of energies, as a function of temperature and as a 
result of a phase transition. A considerable variation is 
expected as a function of electron concentration, and it is 
this aspect on which we have focused in the dimer calcula- 



tions presented here. If the electron concentration is low, 
the pseudo-spin degree of freedom will be rarely excited, 
and the electrons will be highly coherent. However, if the 
electron concentration is high, then electron movement 
will be accompanied by pseudo-spin excitations. There 
is considerable experimental evidence for such incoherent 
processes in the cuprates, namely the mid-infrared band. 
Our calculations clearly indicate that the Drude-like por- 
tion for holes has reduced mobility compared to that of 
electrons. The connection of the optical sum rule to the 
kinetic energy, and how this probe can demonstrate this 
physics has been worked out in great detail for the dimer 
system considered here. More detailed comparison to ex- 
periment will have to rely on DMFT calculations^ that 
provide answers in the thermodynamic limit. 

The results of the dimer calculations presented here 
agree with the physics originally obtained in a model in 
which the pseudo-spin degree of freedom impacted the 
on-site energy of an electron.^ Here, the pseudo-spin de- 
gree of freedom alters the effective Coulomb interaction 
between two electrons through a dynamical change in 
the on-site electron-electron interaction epitomized by 
t/E5l The qualitative picture obtained here also provides 
a better understanding of the conclusions obtained for 
an infinite lattice in Ref. |16j : holes are less mobile than 
electrons, and the optical spectral weight distribution is 
significantly different for holes than for electrons. 
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